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Abstract. We present a systematic study of elliptic flow as a function of centrality, 
pseudorapidity, transverse momentum and energy for Cu+Cu and Au+Au collisions 
from the PHOBOS experiment. New data on elliptic flow in Cu+Cu collisions at ^/s NN 
= 22.4 GeV are shown. Elliptic flow scaled by participant eccentricity is found to be 
similar for both systems when collisions with the same number of participants or the 
same average area density are compared. This similarity is observed over a wide range 
in pseudorapidity and transverse momentum, indicating that participant eccentricity is 
the relevant quantity for generating the azimuthal asymmetry leading to the observed 
elliptic flow. 
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1. Introduction 

The characterization of elliptic flow has proven to be one of the most fruitful probes 
of the dynamics of heavy ion collisions at RHIC It originates from the almond shape of 
the overlap zone of the collision which produces, through unequal pressure gradients, an 
anisotropy in the transverse momentum distribution [1]. The dominant contribution to 
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Figure 1. Panel a) elliptic flow (V2) vs. r\ for Cu+Cu collisions at ^/s NN — 22.4 
GeV for centrality 0-40% using the hit-based method. The boxes show the systematic 
errors. Panels b), c), d) and e) V2 vs. p T for Au+Au and Cu+Cu collisions at ^/s NN 
= 62.4 and 200 GeV, b) and d) for centrality 0-20% , c) and e) for centrality 20-40%. 
The bars in the plots represent the statistical errors. 



this anisotropy is due to elliptic flow and is measured by the second coefficient, v 2 , of the 
Fourier expansion of the azimuthal distribution of produced particles. The large value 
of v 2 observed experimentally in semi-central Au+Au collisions at RHIC is consistent 
with non- viscous hydrodynamic expansion of quark gluon plasma (QGP) droplets [2]. A 
strong pseudorapidity dependence of elliptic flow reported by PHOBOS [3-5] provides 
useful information for constraining models of the full three-dimensional hydrodynamic 
evolution of the system. 

In this paper, we present elliptic flow of charged hadrons in Cu+Cu and Au+Au 
collisions at v /s NN = 19.6, 22.4, 62.4 and 200 GeV as a function of pseudorapidity, 
centrality and transverse momentum. The measurements of elliptic flow in 22.4 GeV 
Cu+Cu collisions are shown for the first time. This work completes our systematic study 
of elliptic flow measurements, providing an extensive and precise set of experimental data 
for Cu+Cu and Au+Au collisions at RHIC. Furthermore, the comparison of the data 
from Cu+Cu and Au+Au collisions measured by PHOBOS experiment provides new 
information on the interplay between initial state collision geometry and elliptic flow. 

2. Results and Initial Eccentricity 

The Cu+Cu and Au+Au data presented in this work were analyzed in the same 
way, using the "hit-based" and "track-based" analysis methods [3]. Fig. [Th, shows 
the preliminary results of the elliptic flow signal as a function of pseudorapidity (rj) 
in the Cu+Cu collisions at y/s NN = 22 .4 GeV for 0-40% most central events. The 
Cu+Cu v 2 displays a strikingly similar shape in rj to Au+Au collisions at nearly the 
same energy (19.6 GeV) [6]. The strength of Cu+Cu v 2 signal is surprising in light of 
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Figure 2. e par t and e s td of the collision zone of Cu+Cu and Au+Au at ^/s NN =200 GeV 
as function of N part from PHOBOS Glauber MC. The continuous and dashed lines 

correspond to the RMS of £ par t(\/ (£ par t))- The gray bands are discussed in the text. 
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Figure 3. Panel a) shows V2/(e par t) as function of mid-rapidity (|^| < 1) particle 
area density l/(S)(dN/dy) for Cu+Cu and Au+Au collisions. Panel b) v 2 /(£ par t) as 
a function of p T for Cu+Cu and Au+Au collisions at 200 GeV with the same area 
density (same (N part ) = 82). The bars in the plots represent the statistical errors. 



expectations that the smaller system size would result in a much smaller flow signal [7]. 
The dependence of v 2 on the transverse momentum (p T ) of charged hadrons in Au+Au 
and Cu+Cu collisions at v /s^7 = 62.4 and 200 GeV for centrality bins 0-20% and 20-40% 
are presented on Figs. [Tb, lc, Id and le. We observe that for both collision systems, the 
dependence of V2 on p T is similar for the two measured centrality classes. For a given 
system, higher values of v 2 (p T ) are observed for more peripheral collisions. 

In order to distinguish collision dynamics from purely geometrical effects, it has 
been suggested that the measured v 2 should be scaled by the eccentricity of the nuclear 

2 2 

overlap area [8]. The standard definition of the eccentricity is, e st d = I, 1 , where 
o~1 (<Ty) are the variance of the participant nucleon distribution projected on the x (y) 
axis, taken to be along (perpendicular to) the impact parameter direction. 

It has been shown that the measured v 2 in Cu+Cu collisions at RHIC [4, 5] is 
surprisingly large even for most central collisions, for which the average eccentricity 
of the overlap region is small. The PHOBOS collaboration has shown that for small 
systems or small transverse overlap regions, event-by-event fluctuations in the shape of 
the initial collision region affect the elliptic flow. Monte Carlo (MC) Glauber studies 
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Figure 4. Panels a) and b) show the v 2 /(£ P art) as function of 77 for Au+Au (Cu+Cu) 
collisions with same (N part ) at 200 and 62.4 GeV, respectively. Panel c) shows 
v 2/(£part) as a function of rj — \r\\ — ybcam for Cu+Cu and Au+Au collisions. The 
bars in the plots represent the statistical errors. 



have shown that the fluctuations in the nucleon positions frequently create a situation 
where the minor axis of the overlap ellipse of the participant nucleons is not aligned with 
the impact parameter vector. To account for this effect, PHOBOS has introduced the 

participant eccentricity defined as [4,5]: £ part = — — y 2 t 2 — — , where a xy = (xy) — (x)(y) 
is the covariance. This definition accounts for the nucleon fluctuations by quantifying 
the eccentricity event-by-event with respect to the overlap region of the participants 
nucleons. Fig. [2] shows the Glauber model calculations of e st d, £ pa ,n and J (£p art ) as a 
function of N part for Cu+Cu and Au+Au collisions at =200 GeV. The gray bands 
correspond to systematic errors obtained by varying the Glauber model parameters 
such as the nuclear radius, nuclear skin depth, nucleon-nucleon inelastic cross-section 
and minimum nucleon separation. We observe that e part and y (£p ar t) distributions 
are similar, within the small systematic errors, for both systems. The importance of 
the definition of eccentricity in comparing Cu+Cu and Au+Au results is presented 
on Figs. [3] and 4, showing the eccentricity scaled elliptic flow, v 2 / (e par t), for the two 
collision systems. For the comparison we selected centrality bins in Cu+Cu and Au+Au 
such that (Np art ) are matched. For such centrality bins also the average area density, 
1/(S) (dN/dy), is approximately the same. We observe in Fig. [3k, that the v 2 scaled 
by £ P art are similar for both Cu+Cu and Au+Au collisions at the same value of the 
average area density (similar (N part )). It should be noted that in Fig. [3^ which has been 
introduced previousely in Ref. [9], in the y-axis the V2(ry) has been converted to v 2 (|/) 
by scaling the data by factor 0.9 and also in the x-axis the dN/dy = 1.15 dN/d?7 at 
mid-rapidity region, \r)\ < 1. This similarity between Cu+Cu and Au+Au collisions in 
Fig. [3k. is also observed as a function of transverse momentum (see Fig. [3b) as well as 
in a wide pseudorapidity range as shown in Figs UK and 4b. Furthermore, Fig. |4b shows 
that the Cu+Cu and Au+Au systems at 62.4 and 200 GeV exhibit the same extended 
longitudinal scaling when e part and (N part ) are taken into consideration. It should be 
noted that within experimental errors, similar scaling properties should be observed 
using a (£part); as advocated in Ref. [10]. 
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3. Summary 

In summary, we have performed a comprehensive examination of the elliptic flow 
of charged hadrons produced in Cu+Cu and Au+Au collisions at y/s NN = 19.6, 22.4, 
62.4 and 200 GeV as a function of pseudorapidity, centrality and transverse momentum. 
The measurements of elliptic flow in 22.4 GeV Cu+Cu collisions are shown for the 
first time. The comparison of the data from Cu+Cu and Au+Au collisions provides 
new information illustrating that the participant eccentricity is the relevant geometric 
quantity for generating the azimuthal asymmetry leading to the observed elliptic flow. 
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